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Abstract. Various alternative theories of gravity predict dipolar gravitational 
radiation in addition to quadrupolar radiation. We show that gravitational wave 
(GW) observations of inspiralling compact binaries can put interesting constraints on 
the strengths of the dipole modes of GW polarizations. We put forward a physically 
O . motivated gravitational waveform for dipole modes, in the Fourier domain, in terms 

of two parameters: one which captures the relative amplitude of the dipole mode with 



^ ■ respect to the quadrupole mode (a) and the other a dipole term in the phase (/3). We 

then use this two parameter representation to discuss typical bounds on their values 
using GW measurements. We obtain the expected bounds on the amplitude parameter 
a and the phase parameter /3 for Advanced LIGO (AdvLIGO) and Einstein Telescope 
(ET) noise power spectral densities using Fisher information matrix. AdvLIGO and 
ET may at best bound a to an accuracy of ^ 10^'^ and ^ 10^^ and /? to an accuracy 
of ~ 10~^ and ~ 10~^ respectively. 
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1. Introduction and summary 

General Relativity (GR) is the most empirically successful theory of relativistic gravity 
to date. Despite its success, there is strong motivation to test it in regimes where it 
has not been tested. Every relativistic theory of gravity has its own predictions about 
the nature and properties of gravitational waves (GWs). Hence observations of GWs 
from various astrophysical phenomena can be used to test various theories of gravity 
including GR (see e.g. [H Ej |3]). Among the various sources of GWs, inspiralling 
compact binaries are the most promising, both for the first detection as well as for 
testing alternative theories of gravity. This is because, for these sources one can very 
accurately compute the gravitational waveforms within GR. 

From a GW point of view, the modifications of GR can happen in two ways: one 
is in the generation of the waves, because of the presence of other fields like scalars or 
vectors which couple to the metric; and the second is in the effects of propagation 
of the waves. The major difficulty in analysing the effects of alternative theories 
of gravity in the context of GWs is that unlike GR, where very precise waveforms 
are available for inspiralling compact binaries (see e.g.jU |5]), theoretical progress in 
modelling the GW emission within the framework of alternative theories of gravity 
is scanty, apart from scalar-tensor theories (see e.g. [6] for a review of Chern-Simons 
theory and [3 E] for some recent works on scalar tensor theories). Hence one has to 
come up with parametrized descriptions of the waveform within different alternative 
theories of gravity. For example. Will discussed model gravitational waveforms in the 
case of scalar tensor theories [9] (see also [10] ) and massive graviton theories [11] . The 
waveforms could be expressed in terms of a one-parameter deviation from GR in the GW 
phasing in both cases. He also discussed the expected bounds on these parameters from 
various GW observations. These estimates were revised later accounting for various 
physical effects such as spin and using more up-to-date waveforms and detector noise 
characteristics [T0l[I3[T3l[Tll[T5l[T6l[I71[T8l[l9]. 

There have also been attempts to write down generic gravitational waveforms 
in terms of the differences in the phasing coefficients, and to study the ability of 
GW detectors to measure a few of them. This amounts to performing consistency 
tests of post-Newtonian (PN) theory within GR in the mass plane of the binary 
components [201 IHl |22]. More recently, Yunes and Pretorius proposed a general 
parametrization of gravitational waveforms called the Parameterized-Post-Einsteinian 
(ppE) formalism [231 |2l]. 

A generic metric theory of gravity has six modes of GW polarizations as opposed 
to the two predicted by GR [251 EE]. If these additional modes were to be detected 
by GW experiments, it would rule out GR. Alternatively, GW observations can put 
stringent constraints on the parameters governing these additional polarization modes. 
If no assumptions are made about the source geometry or the emitted waveform, one 
requires six independent GW data streams plus knowledge of the source direction to 
measure or constrain all six modes of GW polarizations. This would mean a network 



Generic hounds on Dipolar gravitational radiation 



3 



of six independent ground based-GW detectors which does not seem plausible in the 
immediate future. 

In this paper we show that one can constrain a restricted class of modes 
of polarization, which correspond to dipolar GW emission, using interferometric 
measurement of the GW signal. We discuss the bounds on a class of theories which 
predict dipolar GW radiation in addition to the quadrupolar GW emission similar 
to GR. We propose a physically motivated waveform model of dipolar GW modes of 
polarization, containing two free parameters, that could arise in a theory of gravity 
with a vector field or scalar field in addition to the usual metric. We then use this 
two-parameter description (in addition to the standard GR parameters) to discuss the 
bounds possible from ground-based GW interferometers. In addition to scalar-tensor 
theories like Brans-Dicke theory which predicts dipolar GW emission, the recently 
discussed theories with vector fields like the TeVeS theory [271 EH] and Einstein- Aether 
theory [29], and the Modified Gravity theory (MOG) [30] also predict dipolar GWs. 

The paper is organized in the following way. Sec. |2] lists the assumptions we 
are making to obtain a gravitational waveform, including the contributions from the 
dipolar modes. The Fourier- domain waveform is derived in Sec. [3l The details of the 
Fisher matrix analysis are presented in Sec. HJ and the results are discussed in Sec. |5l 
Conclusions are provided in Sec. |6l 

2. Assumptions 

In order to write down a parametrized waveform which captures both the quadrupolar 
and dipolar GWs, we make the following set of assumptions. 

(i) For the quadrupole part we use the GR waveform. In other words we completely 
ignore the corrections to the quadrupolar component of the waveform from the 
underlying theory of gravity. This is consistent with the theme of the paper where 
we focus on the dipole component of the gravitational waveform. One can include 
the corrections to the GR waveform following the general parametrizations given 
in [20l [23] which is outside the scope of this paper. 

(ii) The phasing of the dipolar GW mode will be proportional to the orbital phase 
instead of twice the orbital phase as for the quadrupolar tensor mode. Further, 
because of a dipole component in the energy flux, the orbital phase evolution will 
also have dipole term similar to the case of scalar-tensor theories [9]. However, we 
will ignore higher harmonics of the orbital phase above the quadrupolar one. 

(iii) From dimensional arguments, the frequency dependence of the phasing term 
from the dipole contributions should be 0(w~^) relative to the GR quadrupole 
contribution, and that of the amplitude should be 0{ir^) times the leading 
quadrupole term, where v represents the characteristic source velocity. This is 
easy to understand intuitively, as in the total GW flux for the quadrupole and 
dipole terms differ by {d/dtY ~ 0{v'^) whereas for the waveform they differ by 
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d/dt ~ 0{v). (For the multipolar structure of energy flux and waveform see Eq. 
(12.2) of Ref. [31] and (3.11) of [32] respectively. Though these papers discuss 
multipoles higher than the quadrupole, it is easy to infer the v dependence of the 
dipole term from them.) 

(iv) In the underlying theory which predicts dipolar radiation, there could also be 
monopole radiation. However, by analogy with scalar-tensor theories, it is likely 
that such monopole contributions will begin at post-Newtonian order, since the 
monopole (the mass) is conserved at Newtonian order, and will thereby have the 
same scaling with v as the quadrupole contributions. We have also not considered 
possible deviation of effective stress energy tensor of GWs from that of GR as 
discussed in Ref. [33] . 

(v) We have assumed that the dipolar modes propagate at the velocity of light. As 
mentioned earlier, introducing an additional massive graviton parameter [TT], one 
can probe the non-null propagation of these modes. But we avoid this added 
complication in the present work. 

(vi) In the waveform that we write down below, we introduce two parameters a and 
/3, which capture the relative amplitude of the dipolar mode with respect to the 
GR mode and an additional dipole term in the phasing formula, respectively. It is 
likely that either or both of these parameters will depend on the masses and on the 
parameters determining the structure of the compact objects, as in Brans-Dicke 
(BD) theory (see Eq. (90) of [M] for the case of Einstein- Aether theory). However, 
we will not attempt to include such dependences in the parameters a and /3. But it 
should be noted that for objects like neutron stars, these parameters may depend 
upon the equation of state. 

(vii) Lastly, we have considered here single interferometric GW observations. But it can 
be extended to the case of GW detector network. 



3. Derivation of the phasing formula for dipolar modes 

We start with a schematic expression for the detected strain at the detector output after 
projecting onto the response of the detector combining the GR and the dipole modes, 
writing them as 

Kt) = hGK{t) + hY,,^{t), (1) 

= AGR(t)e2^'^W+ADip(t)e'*W, (2) 

= Av^ (t) e2^"^(*) +Aav{t) e'*(*) , (3) 

where \E'(t) is the time domain GW phase, AQii{t) and AY)ip{t) are the time domain 
amplitudes of the GR and dipole modes, a is the parameter which captures the relative 
strength of the dipole mode with respect to the GR mode and v{t) is the characteristic 
instantaneous velocity of the binary. Since the analysis is for a single detector, we have 
averaged over the antenna pattern functions. This averaging yields numerical factors 
which have been absorbed in the deflnition of the two amplitudes. 
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Using the stationary phase approximation (SPA) (see [551 ES] for detailed 
discussions and examples) one can obtain the Fourier-domain gravitational waveform. 
We follow the procedure developed by Van Den Broeck and Sengupta[37] to address the 
complication in computing the Fourier transform when higher multipolar components 
of the binary are to be taken into account. As a simple first step, we include only the 
dipole and quadrupole harmonics, and the SPA leads to the expression 



'2F{V2) \/F{v,) 

where Vk = (27rm //A;)^^^, and F is the frequency sweep, which can be written in terms 
of the GR contribution Fqr plus a dipolar correction. Introducing the parameter /3 for 
the dipolar correction, we express F in the form 

F{vk)=F^r [^ + ^^k' + 0{vl)] , (5) 

where k = 2 for the quadrupole GR mode and /c = 1 for the dipolar modes. Since 
we use the restricted waveform with Newtonian amplitude only the leading term in the 
Fqr expression is required for the calculation, which we denote by Fq™'^. Substituting 
Eq. ([5]) into Eq. (jlj) and keeping terms in the amplitude to linear order in a and (3, we 
obtain 



.(6) 



In the above expression, the Fourier Domain phasing is given by 

^(^'.) = ^GrM(i+^%' +■■■). (7) 

where it is straightforward to show that /3 = —4/3/7. The above expression is normalized 
to the Newtonian contribution to the phasing formula, ^I^giT* ^"^^ account for the 
full GR phasing formula up to 3.5PN [381 [39] in our calculations. 

There is an important structural difference between the two terms which represent 
amplitude deviations in the final waveform. The term proportional to a comes purely 
from the time-domain dipole mode polarization and its relative strength with respect 
to the GR polarization modes. The second term proportional to /3 arises from the 
calculation of the Fourier domain waveform using the SPA (via F) and this parameter 
is related to the strength of dipolar radiation component (energy flux) with respect 
to the quadrupolar radiation of GR. This parameter modifies the phasing of both the 
quadrupole and the dipole modes, as the effective energy balance will involve both the 
dipolar and the quadrupolar components of energy flux. 

It is useful to see how this waveform maps onto the more generalized 
parametrizations mentioned earlier. There are two leading contenders for general 
parametrizations of inspiral waveforms: the 'parametrized tests' of PN theory (PTPN) 
by Arun et al [20| [2T| [22] and the Parametrized post-Einsteinian (ppE) representation 
by Yunes and Pretorius [231 EH- The proposal by Arun et al is to determine, from 
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GW observations, various phasing coefficients and check for their consistencies in the 
mass plane. Their proposal, so far, has not been extended to include the deviations in 
amplitude w.r.t GR. Hence all that we can say is that there will be additional dipole 
contributions to the parametrized phasing expression in their formalism. 

On the other hand the ppE accounts for typical amplitude and phase deviations 
from GR in their parametrization. For this comparison, we rewrite Eq. ([6]) as 



where /3 = \^q^*^ and 6h{vi) is the amplitude term proportional to a in Eq. ([6]). 
Comparing the above equation with Eq. (1) of Ref. [23] and ignoring the merger and 
ringdown part of the waveform, it is straightforward to see that the ffist term of Eq. (Q 
can be written down in a ppE form by choosing OppE = — f?7~^^^, /3ppE = 2/3?7~^/^, 
a = —2/3 and b = —2/3. In the above expression, r] refers to the symmetric mass ratio 
which is the ratio of reduced mass to the total mass. However second term of Eq. ([6]) 
cannot be mapped to any corresponding term in ppE. This is because ppE, in its current 
form, does not incorporate harmonics other than the leading quadrupolar one. In other 
words, the additional terms which arise here should be added to the ppE waveform when 
non-quadrupolar modes are included into their framework. 

4. Fisher matrix analysis 

Our goal now is to obtain the bounds on a and /3 that would be possible with various 
GW detectors and for different types of sources. Using the waveform in Eq. ([6]), we 
perform a Fisher matrix analysis of the proble by treating the following parameters 
of the signal as independent: 



where, as usual, Dl is the luminosity distance to the source, t^, (pc are the time and phase 
of coalescence, and Mc and rj are the chirp mass and symmetric mass ratio of the binary. 
(For details of the parameter estimation using GR waveforms, see e.g. Ref. [H].) The 
errors on a and /3, obtained from the covariance matrix constructed, translate directly 
into the typical bounds on these parameters using various interferometers. For the 
phasing formula, we use the restricted 3.5PN waveform obtained using PN theory [^B] 
and add the dipole term we propose to the phasing. We ignore the effect of spins, 
eccentricity and that of harmonics higher than dipole and quadrupole in this work, as 
a first approximation. 

X See Ref. [40j for one of the first applications of Fisher matrix in GW context and the basic theory 
involved. 





(9) 



r = {«,/?, 0c, log M„ log r/}. 



(10) 




Figure 1. Bounds on a and /? as a function of the mass of the binary system for a 
mass ratio of 2. A lower cut-off in frequency of lOHz and 20Hz has been used for ET 
and AdvLIGO, respectively. 



For the noise spectrum of AdvLIGO detector, we have used the analytic fit given 
in Eq. (2.1) of [22]. The configuration and the noise model for third generation Einstein 
Telescope is based on. Eq. (2.2) and (2.3) of [22] . 

The upper frequency cut-off of signals is assumed to be the frequency at the last 
stable orbit -Flso after which PN approximation ceases to be valid. The quadrupolar 
term is truncated at 2Flso and the dipole term at -Flso; respectively. The expression 
for -Flso as a function of the total mass m of the binary, is given by 

FLso = (2vr63/2^)"^ (11) 

We take the low frequency cut-off of AdvLIGO to be 20Hz and of ET to be lOHz. The 
distance to the source in both cases is fixed to be 200 Mpc. The effect of increasing 
distance will be to increase the errors in a linear fashion. The mass ratio of the systems 
is assumed to be 2. 

One of the caveats of the present analysis is the use of Fisher matrix as an error 
estimator. Fisher matrix may underestimate the errors for low signal to noise ratio 
events [^ HH |15]. There have been different proposals in the literature to go beyond 
the Fisher matrix [46l HT] in the context of parameter estimation of inspiral signals. In 
the present work we do not address these issues and leave them for future. 



5. Results 



Fig. [H presents our results. The left panel shows the results for the bounds on a and 
the right panel for /3 for ground-based interferometers AdvLIGO and ET. The bounds 
obtained from the Fisher matrix analysis are shown as a function of the total mass of 
the binary. The masses range from the binary neutron star or NS-stellar mass black 
hole (BH) case on the lower side to the black hole binaries consisting of stellar mass 
BHs or intermediate mass BHs (IMBH) on the higher end. 
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The Einstein Telescope would be able to constrain the dipolar amplitude to ~ 10~^ 
whereas AdvLIGO would bound this amplitude to be ~ 10"^. The dipole parameter in 
the phasing may be best bounded to a value ~ 10~^ and 10^^ by ET and AdvLIGO, 
respectively. This corresponds to binary neutron star systems or NS and stellar mass BH 
systems. However, it should be borne in mind that the event rates for IMBH systems 
may be far too small at 200 Mpc for AdvLIGO and ET. Hence the best bounds on a 
may be weaker by about 5 times if we put the sources at, say IGpc, where there may 
be more sources. The factor of 5 simply follows from the fact that errors get worse 
inversely with distance. However for the best bounds on /3, which come from binary NS 
systems or NS-stellar mass BH systems, the estimates based on a luminosity distance 
of 200 Mpc may not be unrealistic. 

From the figures its clear that /3 is always much better estimated than the parameter 
a. This can be easily understood because (3 appears in the phasing, which is more 
important than the amplitude as far as parameter estimation is concerned. Though 
there is a /3 dependent contribution in the amplitude also, coming from the re-expansion 
of F, it adds much less to the already dominant effect of (3 in the phase. 

6. Conclusion 

We proposed a method to constrain the amplitudes of the dipole modes of GW 
polarization of a generic metric theory of gravity using parameter estimation on detected 
GW events. In the absence of gravitational waveforms for dipole modes predicted by 
specific theories of gravity, we wrote down a physically motivated, phenomenological 
representation of the waveform. Based on a Fisher-matrix analysis, we conclude that 
Advanced ground-based GW detectors like AdvLIGO and ET would provide interesting 
constraints. 

All these analyses have been performed assuming a single detector configuration. 
Using a network of ground-based GW interferometers, one could improve the estimates 
quoted here. We have not included the effects of higher harmonics, which might again 
improve the estimates. It would also be interesting to study the effect of spins and 
eccentricity on these estimates. We postpone these issues for a future work. 
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